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(— i . Abstract 

A mean-field model of Ising spin glass with the Hamiltonian being a sum of the infinite-range 
ferromagnetic and random antiferromagnetic interactions has been studied. It is shown that this 
model has phase transition in external magnetic field into inergodic spin glass phase with a number of 
metastable states. The thermodynamic properties of metastable states are studied at T = and near 
the transition. The relations between the characteristics of slow nonequilibrium processes in spin glass 
phase (such as hysteresis loop form, thermo-remanent and isothermal remanent magnetizations, field- 
cooled and zero- field-cooled thermodynamic quantities) and thermodynamic parameters of metastable 
states are established. 
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The possibility of existence of spin glass phases as specific thermodynamic states in solid solutions of fer- 
romagnets and antifcrromagnets has been first suggested in pioneering work by Edwards and Anderson 
J]]]. According to|l| , these phases are characterized by the appearance of local spontaneous magnetic 
moments with chaotic orientations, which are determined by a random distribution ferromagnetic and 
antiferromagnetic interactions throughout the crystal. The first attempt to give the quantitative de- 
scription of spin glass transition was made by Sherrington and Kirkpatrick (SK)pfl, who considered a 
mean-field model with infinite range random interaction. But the solution they got appeared to be un- 
Q^ \ stable in the glass phase for T < T sg , H < Hat ~ (T sg — T) ' ||. The further attempts to describe the 

thermodynamics of spin glass phase below the Almeida-Thouless (AT) line, H < Hat, resulted in the 
construction of 'replica symmetry breaking' scheme by Parisi Q being the procedure of the analytical 
continuation in the replica method used in the studies of the SK model. Now it is common belief that 
Parisi solution is stable below AT-line and it is the basic result in the spin glass theory. With some 
reservations concerning the mathematical foundations of the replica method and the suggested procedure 
of the analytical continuation, one may consider the Parisi's solution as the first exact description of the 
O ■ thermodynamics of inergodic spin glass phase. 

Later, the other exactly solvable spin glass models have been studied with the use of the replica 
method (the vector models [pi 6], p-spin spherical models [M 8]) and without it (ordinary spherical spin 
glass [pi, Bethe lattice spin glass pi). Let us also notice the study of the mean-field equations for the local 
magnetic moments of the SK model pd[ |. In these studies the equilibrium thermodynamic parameters 
averaged over random interactions were obtained in the framework of the standard statistical mechanics. 
Still these results appear to be insufficient for the description of properties of real spin glasses. The 
reason lies in the known inergodicity of the spin glass phases Jfj|, i. e. the existence a large number 
of metastable states in these phases. So the experiments give the physical quantities proper to one of 
these states, in which crystal comes depending on a regime and a sequence of cooling and application 
of magnetic field |13], 17]. Meanwhile, the equilibrium quantities refer only to the state with the lowest 
thermodynamic potential and could be obtained after sufficiently large observation time such that crystal 
could come to the lowest equilibrium state when field and/or temperature are changed. As the barriers 
between metastable states are macroscopic (divergent as N — ► oo), the corresponding relaxation times 
are generally astronomically large. Thus crystal can stay in the initial state while his potential becomes 
larger then some other state. 

The situation can be elucidated by considering the uniaxial (Ising) ferromagnet being the simplest 
inergodic system below T c . In the fields, smaller then coercive one, it has two stable states: equilibrium 
state with magnetization m parallel to H, and metastable one having m antiparallel H and greater 



potential. The standard result of the equilibrium statistical mechanics for the dependency m on H in 
this case is a function with a jump at H = corresponding to the equilibrium states, while in the 
real experiments the hysteresis loop is observed, in which metastable states are also present. Thus the 
description of real experiments needs also the description of the properties of metastable states. In this 
example this is trivial at least in the mean-field approximation, but in spin glasses the description of the 
properties of a large number of metastable states is rather difficult task. In particular, the description of 
magnetic properties of a spin glass phase must include a set of functions m(H), corresponding to various 
metastable states, and their lines on the m — H plane would fill some region around the origin. In a 
rather simplified form, the theoretical problem is to determine the boundaries of this inergodic region 
and the values of various thermodynamic quantities for all points inside it. Such theoretical results could 
describe a number of slow nonequilibrium processes in the inergodic spin glass phases. 

Meanwhile, all data about metastable states obtained via Parisi's ansatz come to the probability 
distribution of the 'overlaps' of magnetizations in various metastable states [fill. Now it is not clear 
how this distribution can describe the real experiments. Also there was not got any information on the 
properties of metastable states of the SK model from the TAP equations [0|, it was only established 
that their number is exponentially large |H . 

The studies of the other spin glass models have not also result in the elucidation of the physical 
characteristics of metastable states. The only exception is the Ising spin glass on the Bethe lattice for 
which the numerical study of the internal field distribution has explicitly shown the existence of a number 
of metastable states at T — ||10f| - 

Generally speaking, the study of the properties of metastable states is not necessary for a description 
of experiments in inergodic systems as it could be obtained from the study of nonequilibrium dynamics 
which incorporates automatically the effects of their existence. Thus one can get the hysteresis loop as 
a reaction on the large slow varying field in the dynamic treatment of the uniaxial ferromagnet. But in 
doing this one must eliminate large unobservable times of the relaxation between metastable states. Such 
elimination procedure has been developed in the study of Langevin dynamics of SK model in zero field 
Q , which help to establish the difference between the unobservable equilibrium susceptibility and those 
in the local metastable states measured at real times smaller than that of relaxation over macroscopic 
barriers. Still now it is not clear how this procedure could be generalized to describe in the SK model 
the reaction on the alternating field of finite amplitude. Now the inergodic effects in finite fields were 
described only in the simplest spherical spin glass model having two stable states |li| • 

To resume we may say that the spin glass models considered until now are not sufficiently simple to 
obtain the thermodynamic properties of metastable states needed for description of inevitably nonequi- 
librium processes in spin glasses. Hence it seems worthwhile to find and to study more simple models 
allowing for more complete description of physical properties in the inergodic glass phases. Here we 
present one such model which allows the analytical description of thermodynamic characteristics of all 
metastable states at T — and near the transition point. 

2 Hamiltonian of the model and its properties at T = 0. 

The most simple treatment of metastable states is possible in the framework of mean-field models with 
infinite range interactions when it comes to the determination of local minimums thermodynamic poten- 
tial. Also the random spin glass interaction should be sufficiently simple in all its realizations to make 
possible an analytical treatment. In Ising spin glasses the glass transition results from the competition 
between ferromagnetic and antiferromagnetic interactions so the task is to imitate this competition using 
simple infinite range interactions. 

Let us consider the system having N Ising spins Si a = ±1, divided on iVj, blocks which consist of N s 
spins, N = NbN s . Here index a is the block number and index i is the number of the spin inside the 
block. So the magnetization of the block is 
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and the total magnetization is 
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Let us also introduce the 'antifcrromagnetic' order parameter 

Ha = m a — to 



(3) 



The model Hamiltonian is a sum of ferromagnetic and 'antiferromagnetic' interactions and the external 
field term 
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Here H is the homogeneous external field, J > - ferromagnetic exchange integral with a fixed value, 
Ji > 0- random 'antiferromagnetic' exchange integral. 

The term proportional to J\ imitates the antiferromagnetic bonds distributed throughout a crystal so 
we assume that A/& >> 1. Also we assume that N s diverges in thermodynamic limit (N — > oo) in order 
the mean-field approximation to be valid. 

In general, for large Ni, the 'antiferromagnetic' term in Eq.H gives rise to a number of various types of 
ordering with m = 0, so it would be more correct to call it the 'glass' term. The more so that the block 
analog of the Edwards- Anderson order parameter 
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can be represented as 
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Thus we may say that for TVf, >> 1 the model Hamiltonian, Eq.|4[ describes the competition of ferromag- 
netic and 'glass' ordering. 

For N s — > oo it is easy to find the (nonequilibrium) thermodynamic potential depending on m a : 

F(m) = H (m) /N - TS (m) 



where S (m) is the entropy per spin: 

5(m) =7V _1 ln 
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For N s — > oo the description of the equilibrium thermodynamics comes to the finding of the lowest 
minimum of -F (m) while the less deep minimums correspond to metastable states. The equations defining 
the extrema of F (m) are: 



Tare tanh m a + (Jj — J) to — 3\m a = H 
The minimums correspond to the solutions of Eqs.R with the positively defined Hessian 



d 2 F (m) 



= [T/ (1 - ml) - Ji] S af3 + (Ji - J) /N 



dm a dm/3 

It is easy to show that this simple model is inergodic at T = 0. In this case Eqs.0 become 

m a — sign [H + (J — J{) m + Jim a ] 



(7) 
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When \H\ > max (J, 2Ji — J) Eqs.^have unique solution m a — signH, while at \H\ > max (J, 2J\ — J) 
they also have a number of solutions with arbitrary signs of m a limited only by the condition 



[H+(J- Ji)to] 2 < J 2 



(10) 
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Figure 1: Field dependencies of magnetization in metastable states (dashed lines) and stable states (solid 
lines) at T=0. 



All these solutions are stable thus corresponding to the metastable states of the model, 
magnetization in these states acquires a set of discrete values 



The total 
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Here n is the number of blocks with magnetization ra a = 1, n = 0, ..., Nb. There are ( ra b ) states with a 
given m which differs by the permutations of m a . The total number of inhomogcncous metastable states 
cold be up to 2 Nb — 2 for a given H . 

The energy per spin in these states is determined by their magnetization 



E=(J - J 1 )m 2 /2 - mH - J x /2 

and the entropy, Eq.|6| is zero. 

The equilibrium magnetization corresponding to the states with minimal energy is m eq (H) 
for J > J\ and 
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when J < J\. Here e n = j H _j — j^- + 1, and 9 is the Haviside's step function. The field dependencies 
of magnetization in the equilibrium and metastable states are shown in Fig. 1 . The steps of the function 



f eq 



(H) for J < J i demonstrate the existence of the first order phase transitions at field values 



H,. 
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n = 0, ..., TV;,, at which the upturns of the block magnetizations take place. Qualitatively, just this behavior 
of m eq (H) one may expect in the glass phase of the magnet with equal concentrations of ferromagnetic 
and antiferromagnetic bonds, while the case J > J\ gives a picture proper for a magnet with domination 
of ferromagnetic bonds. Further we will consider the most interesting case when model Hamiltonian 
imitates a spin glass, so we assume that probability distribution for J\ 1 P{J\) , is zero for J > J\ and 
has the form 



P( Ji) = 9 ( Ji - J) W (Ji - J) 



(12) 



In general, one may notice that averaging over random interaction is superficial when all thermodynamic 
parameters can be obtained for every random realization as experimental data are not usually averaged 
over a number of samples. Still in more complex models it is often possible to find only average equilibrium 
quantities. So it is interesting to compare them with the corresponding quantities in the metastable 
states. Thus in present model the magnetization in all states does not depend on the field and magnetic 
susceptibility in them is zero at T — 0. In the same time, averaging of Eq.[fl] over Ji gives for Nf, >> 1 
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As one can see from Fig. 1(a), the non-zero (Xeq) appears due to the series of transitions between 
metastable states. Indeed the differentiation of EqJO gives 



Xeq = 



dm eq (H) 
dH 
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When Nb — ► oo EqJT^ becomes in the sense of distributions 

1 



Xeq 



(J, -J) 



(Ji - jy 



H' 



(15) 
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The average value of this expression coincides with Eq.[f4|. Thus (xeq) is generally the unobservable 
quantity as it describes the changes of m in the series of transitions which take place only astronomically 
large time scale of overbarrier relaxation. What is more, there are no traces of the transitions in the 
limiting expression in Eq.[l6| when the number of metastable states goes to infinity (Nb — > oo) in the 
thermodynamic limit. The indication of their presence via delta functions in \eq exists only when this 
number stays finite in the limit N — > oo. 

Meanwhile in the framework of the present model \eq from EqJIq contains some information on the 
reaction of inergodic system on the slow varying external field. Thus the application of slow AC field with 
amplitude greater than 2 J\ — J would give a hysteresis loop and Xeq defines its slope. In the same time 
ordinary measurement of susceptibility in small fields would give zero value. Possibly the qualitatively 
similar meaning Xeq has in the real systems. Yet we must note that Xeq in present model is a non-self- 
averaging quantity in the sense that being constant before averaging it becomes a function of H after 
it, see EqJll. Still more dramatic effect is caused by the averaging on the nonlinear susceptibilities. It 
follows from Eq.|l| (k > 1) 
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The absence of self-averaging of magnetic susceptibilities seems to be the specific property of the present 
model in which small fluctuations of J\ can cause large deviations in m eq and could be absent in more 
realistic spin glass models. 

Let us note that singularities of non-averaged susceptibilities at H = ± ( J\ — J) correspond to the 
points of the transitions from the inhomogeneous phase into the homogeneous one, thus J\ — J has the 
meaning of the (non-averaged) Almeida-Thouless field. The corresponding anomalies of the averaged 
quantities would exist at finite H — ±Hat if function W (J') in Eq.[l2| has a bounded support, i.e. when 

W ( J') = for J' greater some J > and W ( J') > otherwise. Then Hat = J and anomalies of (xeq ) 

for H — > ±Hat will be determined by the behavior of W (J 1 ) at J' — > J. This behavior determines also 
how the block Edwards- Anderson order parameter, Eq.||, goes to zero when H — > ±Hat- For Nb — > oo 
we get from Eq]5j 
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(Qeq) = 1 - (m 2 eq ) 
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In general case (W (j) < oo) the transition into spin glass phase with (q eq ) ^ at \H\ < Hat is 
not accompanied with divergencies of ( Xeq ) in contrast with SK model. This is because the upturn 

of the last block along the field is the first order transition. So the instabilities which could cause such 
divergencies does not occur in the present model. 

Let us also note that there are two ferromagnetic phases: inergodic phase with a number of inho- 
mogeneous metastable states and ergodic one having the unique ferromagnetic state. The average value 
of the field corresponding the transition point between these phases can be obtained by considering the 



boundaries (upper and lower) of the region where metastable states exist which represent also the upper 
and lower branches of hysteresis loop . They are (see Fig. 1(a)): 



m±(H) 



H±Jt 



J i - J 
The averaging of this expression over J\ gives 



J 1 -jf-(H±J 1 f + sign(H±J 1 )9 (H ± Jif - ( J x - J)' 



(m±(H)) = sign(H±J) + e(TH-J) f dJ'W {J') ( ^-^- ± 2 J 



(18) 
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For \H\ greater some H e these branches coincide thus indicating the transition into ergodic phase. 
The condition defining H e is vanishing of the integral in Eq.[18| so H e will be finite when W (J') has a 
bounded support. In this case 

H e = 2J + J = 2H AT + J 



There exists a functional relation between (m± (H)) and (m eq (H)), Eq.[13]of the following form: 

(th- j n 
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From Eq.ffCfl follows also: 



d -^^ = e(TH-j)( x J< [f - J 



(19) 
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These relations are specific for the model under consideration but, probably, some relations between field 
dependency of the average equilibrium magnetization and hysteresis loop contour exist also in other spin 
glass models at T = 0. 

For the simple 'rectangular' function W 



we get 
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The divergency of magnetic susceptibilities at H = is a consequence of W (0) ^ 0. They would be finite 
at zero field if W goes to zero as some power of J' or faster when J' — ► 0. For example, the averaging 
with the 'triangle' function 



W (f) = 2J'9 (f - J) /J 2 



(21) 



gives the following results 
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The averaging of the thermodynamic parameters of metastable states existing inside hysteresis loop is 
trivial their magnetizations do not depend on H and susceptibilities are zero. Also q — 1 — m 2 and S = 
in all states. But we must note that equilibrium entropy is not strictly zero as for a given m there are 
{n (l— ml/2/ states with equal potentials so the configurational entropy term 

^-n/ (m)=iV- 1 In (^ {l Nb m)/2/ 
is added to the expression in Eq.g. But in the thermodynamic limit S con f (m) goes to zero as 2V" 1 . 

3 Thermodynamics near the transition. 

The stable inhomogeneous solutions of equations of state, Eq.[7], appear at T < J\. So in the case 
considered here ( J\ > J) there is a second order phase transition from homogeneous paramagnetic phase 
into inhomogeneous inergodic spin glass one at T = J\,H — 0. Let us consider the thermodynamics of 
the model in the vicinity of this transitions assuming 

m Q « 1 (22) 

In this case Eq.ffl acquire the form: 

niria + (t — Ti)m + m a /3 = h (23) 

Here h = H/T, n = 1 — Ji/T, t — 1 — J/T, T\ > r. Hessian, Eq.0, becomes in this region 

T - 1 £3£- (n+<) ^ +{T - n )/N (24) 

It follows from Eq.^|| and Eq.J23| that 

h « l,n << l,r << 1 

For these conditions to be fulfilled for every random J\ we must assume that W ( J') in Eq.[12| has 
sufficiently narrow bounded support, that is the possible values of J\ — J must be less than some J > 
obeying the condition 

J « J 

When Ti > Eqs.g3]have unique homogeneous solution. Let us denote it as mo. It does not depend on 
T\ and obeys the equation 

titio + ml/3 = h (25) 

When n < Eqs|2J| have 3 Nb ~ 3 inhomogeneous solutions beside m^. As all m a obey the same 
cubic equation they can acquire only three different values which we denote as fh s ,s = 1,2,3. Then all 
inhomogeneous solutions can be characterized by the three numbers n s ^ Nb, 



Y,n s = N b 



which show how many m a have the value m s . There are — , — -, — r solutions which differ by the permuta- 
tions of m a and the total number of solutions is 



y- N b 
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But only 2 Nb — 2 of them could be stable. Indeed, Hessian, Eq.g4[ has three eigenvalues equal to t\ + m 2 
with degeneracy n s — 1, which correspond to the eigenvectors having a zero sum of components. There 
are also three non-degenerate eigenvalues which are the solutions of the equation 



1 + ^E +~ 2 A =0 (26) 

Nb '-~ l T\ + mi — A 



Using the Viet's theorem for Eq]2^ according to which 

^m s = (27) 

S 

y rh s rh s > — 3n (28) 

s<s' 

we can get the relation 

*— ' Ti + TO« 

It shows that all three eigenvalues t\ + m 2 . could not be positive simultaneously so the stable solutions 
must have at least one of the numbers n s equal to or 1. But if all n s > then one of the solutions of 
Eq.p6| becomes negative for large Nj, » 1. Thus the stable solutions must have one of n a equal to zero. 
Further we will consider just these solutions putting 713 = 0. The stability condition for them reduces to 
one inequality 

ti + m\ < (29) 

It follows from Eq.B7^ and Eq.Eq that fh s can be represented in the following form 

fhi = 2 (-ri) 1/2 cos Up - -J 
m 2 = -2 (-ri) 1/2 cos Up + -J 

m 3 = -2(-ri) 1/2 sin(^ (30) 

so the stability condition, Eq]29, is equivalent to the inequality 

M<tt/6 (31) 



It follows from the definition of m: 
Here 



m = v\m\ + V2rrii 



n s 
v* = — 

N b 

V\+V2 = ^ 



(32) 



so 



m — (—ti) ' ( V3 A cos <p + sin tp J (33) 

A = v x — i> 2 
Inserting EqJ3^ into Eq.g3| we get 

3(r-ri)m = 3/i-2(-Ti) 3/2 sin3(p (34) 



Excluding m from Eq.p3| and Eq. 34 we obtain the equation for tp 



2 (— n) sin3(/P + 3 (r — n) (— n) ( v3 A cos tp + sin tp ) =3h (35) 



At all t > Ti the left side of Eq]35J is a monotonously growing function of tp for \tp\ < tt/6. Hence, there 
is only one stable solution for fh s at a given A. There are ( n 6 ) metastable states corresponding to this 
solution which differ by m a permutations. 

The explicit solution of Eqj3| can be found for A = when it becomes cubic. In the limiting cases 
A = ±1 Eq.pq reduces also to a cubic one for fh\ or mi which coincides with EqJ2q. 




Figure 2: Field dependencies of magnetization in metastable states (dashed lines) and stable states (solid 
lines) near transition, r < 0. 



In general case Eq.pq and Eq.p3 (or Eq.B4J) give a parametric representation of a dependency of the 
homogeneous magnetization in the metastable states with a given A on r, n and h. The parameter ip 
can be excluded from these equations with the result 



(3A 2 + l) 2 r - 3 (1 - 3A 2 ) (1 - A 2 ) n 



8(9A 2 
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+2V3 A „^ 2 A ^ [(3A 2 + 1) n + 4m 2 ] [- (3A 2 + l) n - m 2 ] 1/2 = (3A 2 + lfh (36) 

From the stability condition, EqM], and Eq. 34 it follows that solutions of Eq.pq is stable in the region 



9[(r-n)m-hf < -4rf 



(37) 



which is the band on the m — h plane. The magnetization is a monotonously growing function of h and 
A inside this band so the field dependencies of magnetization can be represented as a set of uncrossing 
lines bounded from above and below by the mo(h) line as shown in Fig. 2. The other thermodynamic 
parameters of metastable states can be obtained by differentiation of thermodynamic potential which 
near transition has the form 



12F/T = 6 (r - 7i) m 2 + 6n J^ i 

3 

Expressed via if these parameters are 
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q = 3ri (A — l) cos ip 

1 + 2 cos 3ip/ I cos if — V3A sin ip 



T — T\ 

S = ln2 + t\ ( 2 + cos 2ip + V3A sin 2ip) /2 
For the heat capacity we get rather more cumbersome expression 

2ti (v^Asin^ — cos^ — 2cos3(/?) + 3 (r — n) (l — A 2 ) cos ip 



C 



2 (cos <p — V3A sin ip) — 4ri cos tp 



In spite of the absence of explicit expression for ip as a function of h, r and T\ the above formulae allow to 
get some notion about the field and temperature dependencies of these quantities. Thus at the boundaries 
of stability region, \<p\ — ±7r/6 or at 



h = (-n) 1/2 (r - 7i) (3A ± 1) /2 ± 2 (~n) 3/2 /3 
q and x _1 has the lowest values 

x" 1 = Ji - J 

q = 9ti (A 2 - 1) /4 



while the entropy and heat capacity are 



S = 1ii2 + 3t 1 (1±A)/2 
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C = 3(1±A)/2-ti/(t-ti) 

It follows from Eq.B9^ that metastable states exist when 
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(41) 
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When h goes to ±h e the more homogeneous states with A — » ±1 stay stable and their magnetization 
tends to mo (±/i e J = ±2 (— n) . However the limiting values of magnetic susceptibility, entropy 

and heat capacity differ from those in homogeneous state: Xo" 1 = -h ( 1 ~ + TO o), ^o = In 2 — m.g/2, 
C = ml/ (r + m§). 

In the middle of the stability band (at <p = or h = A (— 3ti) (t — 7"i)) we get: g = n (A 2 — l), 
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Ji (1 - 2ri) - J, S = In 2 + 3n/2, C 



inhomogencity when /i — * ±1iat, 
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In this case the diminishing of 



(43) 



?7i, x? $ an d C 1 tend to the corresponding values of the homogeneous phase. 

The Almeida-Thouless field Hat, Eq. [13| determines (to the order A^ -1 ) the point of the transition 
into the homogeneous phase. To show this let us find the values A e<J corresponding to the states with 
the lowest potential. Differentiating F , Eq.|3g over A and using Eqs.g3[ £7] and Eq.gg we get 

dF 3 

— — = Tvn-i (mi - m 2 ) /24 

oA 

d 2 F _T(t~ 7i) (mi - m 2 ) 2 

9A 2 ~~ r 

■ +y -Kg, 

1 t-^t T+rr, 



Thus the states with m^ = or, equivalently, <p = (cf. EqJ3fl) have the lowest potential. One can 
see that EqJ3q has solution <p = when A = h/h^r which is possible at ft 2 < /i 2 ^. When ft 2 > h 2 AT 
F (A) has no minima inside the region A 2 < 1 in which it is defined and the minimal values occur at its 
boundaries for A eq = signH. So the transition into homogeneous state takes place at h = ±h,AT- 

As A is a rational number of the form 2n/Nb — 1 (cf. Eq.p3) it can not be exactly equal to h/h,AT at 
all h 2 < h 2 AT . Hence A eq is defined so that A — Ii/Iiat is minimal and can be represented as 
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Inserting this A eq into Eq.|35| we get the corresponding values of tp eq at h 2 < h 2 AT : 
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Inserting A eq and ip eq into the parametric representations of q and m we obtain the equilibrium values 
of these quantities 
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Differentiating TO eiJ over ft we get the equilibrium susceptibility 



Xeq — 



t — 3ti 



(V ;3 »-" 2 



4ti 



iV h -l 



1 



t + m. 



o 



iV 6 (r -ti)(t-3ti) ^ 
ft 



E'to -1 -^) 



,2 ,^1,, 



A/, 



AT 



The equilibrium entropy can be obtained by the differentiation of the equilibrium potential which to the 
e? order is 



F eq = F[A = h/h AT j - TS con f 

where conhgurational entropy S con f is determined by the logarithm of the number of states with the 
same potential F, 

Sc ^ = 7V_lln U(l-X)/2 
As at T — 0, Sconf is of the order A f ~ 1 and can be neglected so 



S eq = ln2+- Tl $ 
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For the equilibrium heat capacity we get 
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The averaging of these expressions over Ji gives at large N b 
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In derivation of these expression we have used the equivalence of the condition h 2 < h AT and the 
inequality m 2 (h) < mp (Hat) = — 3n where m (ft) is a solution of Eq.£5| such as m (h) h > 0. In its 
turn, from Eq.Eq it follows that the last inequality is equivalent to J\ — J > H/niQ. 

Evidently, the average equilibrium parameters transfer continuously into corresponding values of ho- 
mogeneous phase at H > J mo . Thus the average Almeida-Thouless field Hat is defined as a solution of 
the equation Hat — J mo (Hat) or its equivalent 



ml (Hat) = a 



(44) 



where a = J/T - r = (J + J) /T - 1. 

Let us remind that from Eq.E2 it follows that J >> J. The solution of Eq.E3 exists when a > or 
T>T sg = J+J. For T -> T sg , a « J/J we get 



Hat- J (3a) 



1/2 



and for r < 0, -r >> J/ J 



#at w J (-3r) 



1/2 



Let us further consider the average boundaries of inergodic region on. m — h plane, i.e. the branches of 
the average hysteresis loop 
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and ?71q and m^ are the maximal and minimal solutions of Eq.Eq correspondingly which exist at r < 0, 
4t + TOq < 0- When r > or r < 0, 4r + m§ > then ttIq = mo- 

Eq.E5| shows that the branches of the average hysteresis loop coincide when H 2 > H 2 , H e being the 
solution of equation 

ml (H e ) = 4ct 

For T -> T sg , a « J/J we get 

H e sw 2Jcr 1/2 

and for r < 0, — r >> J/ J, i/ e almost coincide with the coercive field for the homogeneous solution 

H e « 2 J (-r) 3/2 /3 

so the hysteresis loops becomes similar to that of ordinary ferromagnet. 

Let us present the explicit expressions for the functions A± in EqM9 for the 'triangular' function W, 



Eq.0: 
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For the same W the average equilibrium parameters at H 2 < H AT are 
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Let us note once more that the average equilibrium parameters are generally unobservable quantities. 
Probably, the experimental values being rather close to them are obtained after cooling in small external 
fields (field-cooled (FC) regime) for T near T sg 13, 17] when barriers between metastable states are 
relatively small and system could relax into the lowest (or close to it) state at a sufficiently slow cooling. 
In zero field cooled (ZFC) regime when field is applied after cooling below T sg in zero field, the observed 
thermodynamical parameters would differ from equilibrium ones as the system would at first be trapped 
in the state with A = and will stay in it if applied field does not exceed h c = (— ri) (t — 4ti/3) /2, 
cf. Eq.|3!j. Thus at h < h c the ZFC parameters are those of A = metastable states. When applied field 
h > h c the system relaxes into the metastable state at the boundary of stability region (on the lower 
branch of hysteresis loop) having some A > which is a solution of Eq.|39]. Inserting this A in Eq.|I(J and 
Eq.^lj we get the values of thermodynamic parameters the observed quantities would relax to in ZFC 
regime at h > h c . 

Similarly, the parameters of metastable states define the other quantities which are determined in the 
slow nonequilibrium processes in the spin glass phase, such as thermo- remanent magnetization, tuthm, 
which remains after FC process and subsequent switching off the field, and isothermal remanent magne- 
tization, miRM, remaining after ZFC process followed by the application for some time (longer than the 
intravalley relaxation time) an external field |13, 17]. Thus mj^M is apparently nonzero only at h > h c 
and an equation defining it can be obtained by putting h — in Eq.pq and inserting in this equation the 
value of A we get from Eq.EM The equation for mjnm can also be obtained from Eq.pffl by putting in it 
h = and A = min (1, h/hAr)- 



4 Conclusions 

The most remarkable feature of the model considered is the possibility to imitate the properties of 
such complex systems as spin glass with the aid of very simple Hamiltonian. It is common belief that 
the existence of a number of metastable states in spin glasses is caused by the frustration of random 
Hamiltonian, that is the absence of unique spin configuration providing the minimal energy |[L3| . In 
the present model the degeneracy of the ground state results from the permutational symmetry of the 
Hamiltonian instead of frustration. Nevertheless there also exists the transition into the inergodic phase 
and its magnetic properties appear to be very similar to those of real spin glasses including a set of 
transitions between metastable states and inclined hysteresis loop resulting from their presence |L3[ 17]. 
One may suppose that more realistic random Hamiltonians can also have some approximate permutational 
symmetry and corresponding quasidegeneracy of ground state more essential, perhaps, than that caused 
by the frustration. 

We may note that present model implies the definite mechanism of the transitions between metastable 
states in a field, namely, the spin-flop transitions caused by the antiferromagnetic interaction between 
macroscopically large spin blocks. It seems rather probable that in some Ising short-range models of spin 
glasses it is possible to distinguish many clusters with the mostly antiferromagnetic interactions at the 
boundaries and relatively weak interactions inside them. Still it is rather evident that there can be also 
a many more other bond configurations in which degenerate spin configurations can exist at some field 
values. Thus it is hard to say to what extent spin-flop transitions are characteristic for real spin glasses 
and if the Hamiltonian, Eq.EL is a reasonable approximation for some random short-range Hamiltonian 
with some specific type of disorder. Still such possibility seems to be rather probable in view of similarity 
of the properties of the model to those of some real disordered magnets [|13| 17]. 

Finally we may state that in spite of the qualitative nature of the model it allows to get some notion 
about the character of theoretical results relevant for the description of real experiments in inergodic 



systems. It shows how these result may look like and how the thermodynamic parameters of metastatic 
states are related to the characteristics of noncquilibrium processes in spin glass phases. 

This work was made under support from Russian Foundation for Basic Researches, Grants N 98-02- 
18069 and 97-02-17878. 
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